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We construct a model for inelastic neutrino- and electron-nucleon scattering cross sections us- 
ing effective leading order parton distribution functions with a new scaling variable Non- 
perturbative effects are well described using the ^„ scaling variable, in combination with multi- 
plicative K factors at low Q^. Our model describes all inelastic charged lepton-nucleon scattering 
(including resonance) data (HERA/NMC/BCDMS/SLAC/JLab) ranging from very high to 
very low and down to the photo-production region. The model describes existing inelastic 
neutrino-nucleon scattering measurements, and is currently used in analyses of neutrino oscilla- 
tion experiments in the few GeV region. 
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Standard PDFs are extracted from global fits to various sets of deep inelastic (DIS) scattering 
data at high energies and high Q^, where non-perturbative QCD effects are negligible. PDF fits are 
performed within the framework of QCD in either LO, NLO or NNLO. 

In order to use the PDFs at low we use a new scaling variable {^w) to construct effective 
LO PDFs that account for the contributions from target mass corrections, non-perturbative QCD 
effects, and higher order QCD terms. 

Our proposed scaling variable, is derived as follows. Using energy momentum conserva- 
tion, the fractional momentum, ^ carried by a quark in a proton target of mass M is 



22(l + yi + (2Mx)Ve^) 

Here M, is the initial quark mass with average initial transverse momentum Pt, and Mf is the 
mass of the final state quark. Assuming M,- = 0, Pj- = we construct following scahng variable 

_ 2x{Q-+Mr+B) 



Q2[l + VI + (2Mx)2/e2] +2Ax' 

where in general Mf = (except for the case of charm-production in neutrino scattering for which 
M/=1.32 GeV). The parameter A is used to account (on average) for the higher order QCD terms 
and dynamic higher twist in the form of an enhanced target mass term (the effects of the proton 
target mass is already taken into account in the denominator of ^w)- The parameter B is used to 
account (on average) for the initial state quark transverse momentum, and also for the effective 
mass of the final state quark originating from multi-gluon emission. A non-zero B also allows us 
to describe data in the photoproduction limit (all the way down to Q^=0). 

In leading order QCD (e.g. GRV98 PDFs), #2,lo for the scattering of electrons and muons on 
proton (or neutron) targets is given by the sum of quark and anti-quark distributions (each weighted 
the square of the quark charges): 

Our proposed effective LO PDFs model includes the following: 

1. The GRV98 LO Parton Distribution Functions (PDFs) are used to describe =^2,lo(^'2^)- 
The minimum value for these PDFs is 0.8 GeV^. 

2. The scaling variable x is replaced with the scaling variable 

3. We multiply all PDFs by vector K factors such that they have the correct form in the low 
photo-production Umit. We use different forms for the sea and valence quarks. 

T^vector I /^2\ Q j^vector / /^2\ ri y~<2 //^2\i 

'^sea \U ) - q2j^(^ ' '^valenceKU ) - - ^D\U )\ q2 ^ 
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where G/j = 1/(1 + 2^/0.71)^ is the proton elastic form factor. This form for the K factor 
for valence quarks is motivated by closure arguments and the Adler sum rule. At low Q^, 
[1 — Gj){Q^)] is approximately Q^/ {Q^ + 0.178). These modifications are included in order 
to describe low data in the photoproduction limit (0.^=0), where ^2^{x, Q^) is related to 
the photoproduction cross section according to 

4. We freeze the evolution of the GRV98 PDFs at a value of = 0.80 GeV^. Below this Q^, 
^2 is given by; 

^^^^ (x, < 0.8) = K''''""-iQ^)^^^/oi^^, = 0.8) 

5. Finally, we fit for the parameters of the modified effective GRV98 LO PDFs (e.g. ^w) to 
inelastic charged lepton scattering data on hydrogen and deuterium targets (SLAC, BCDMS, 
NMC, HI). We obtain an excellent fit with the following initial parameters: A=0.419, 
5=0.223, and Cvi =0.544, C„2=0.431, and C„,„=0.380, with X^/DOF = 1235/1200. Because 
of these additional K factors, we find that the GRV98 PDFs need to be scaled up by a nor- 
malization factor A^= 1.011. 

We now describe the second iteration of the fit. Theoretically, the Ki factors are not required 
to be the same for the u and d valence quarks or the u, d and s sea quarks and antiquarks. In order 
to allow flexibility in our effective LO model, we treat the Ki factors for u and d valence and sea 
quarks separately. As the predictions of our model are in good agreement with photoproduction 
data, and for much of the resonance region, we now proceed to include photo-production data in 
the fit. In order to get additional constraints on the different K factors for up and down quarks 
separately, we use both hydrogen and deuterium data. 

The second iteration includes the additional photo-production and resonance data in the fit,. 



y2 I f^L—Ehad fp. q2 

j^vector (rp-\ — j^vector (/~p-\ ^ 

^2 ' sea-strangev^ ) ~ f-o , (^vector ' '^sea-up\>i ) ~ Q2\nvector ' 

^ Si ~r ^sea— strange ^ ' ^sea- 
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G2 _|_ f-'vector 
^ ..^ 1/1-1 *- I I I - 1 1 I- i] v2d 

'^valence -downK'd )—^ ^ lAU )\ n2 , r-vector 



The best fit is given by A = 0.621 ± 0.009, B = 0.380 ± 0.004, CUf"' = 0.341 ± 0.007, CUf"' = 
0.417 ± 0.024, Cl^'j"' = 0.323 ± 0.051, C^lf^'' = 0.264 ± 0.015, and c^-^'^'^ = 0.217 ± 0.015. 
The sea factors are C-fi'^„^„=0.621, C~p=0-363, and C~r«„g. was set to be the same as 
CZa-down ■ Here, the parameters are in units of {GeV/cf. The fit X^/DOF =2357/1717, and 
= 1.026 ±0.003. The resonance data add to the jndf because the fit only provides a smooth 
average over the resonances. No neutrino data are included in the fit. 
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For high energy neutrino scattering on quarks and antiquarks, the vector and axial contribu- 
tions are the same. At very high Q^, where the quark parton model is vahd, both the vector and 
axial K factors expected to be 1.0. Therefore neutrinos and antineutrino structure functions for an 
iso-scalar target are given by : 

^2^(x,g2) = i,2[§,,^,(^w,e')+§w^,-(§w,e')] . 

There are two major differences between the case of electron/muon inelastic scattering and the 
case of neutrino and antineutrino scattering. In the neutrino case we have one additional structure 
function {x, Q^). In addition, at low there should be a difference between the vector and axial 
Ki factors due a difference in the non-perturbative axial vector contributions. Unlike the vector ^2 
which must go to zero in the 2^ = limit, the axial part of ^2 is non-zero in the 2^ = limit. 

We account for kinematic and dynamic higher twist and higher order QCD effects in ^2 by 
fitting the parameters of the scaling variable and the K factors to low data for ^2^{x,Q^). 
These should also be valid the vector part of ^2 in neutrino scattering. However, the higher order 
QCD effects in ^2 and x,^'i are different. We account for the different scaling violations in ^2 
and x^-i by adding another correction H{x, Q^) to the expression for x^j, 

The differences between neutrinos and charged lepton scattering are accounted for in the fol- 
lowing expressions: 

x.^^{x,Q^) = 2H{x,Q^) [LiKf'^„qii^,,,Q^)-LjKf'^„qj{^^.,Q^)] 

Where / denotes (valence — up), (valence — down), (sea — up), (sea — down), and (sea — strange). 
Detailed expressions are given in reference [|^]. 

With the above assumptions we calculate the differential cross sections for neutrinos and an- 
tineutrino scattering. We also correct for nuclear effects in iron using the ratio of iron to deuterium 
structure functions as measured in muon and electron scattering experiments. 

Our predictions are in good agreement with the CCFR CDHSW neutrino and antineutrino 
differential cross sections. 
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